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ABSTRACT 



The study of spherical-deformed ground-state phase transitions in fi- 
nite nuclei as a function of N and Z is hindered by the discrete values of 
the nucleon number. A resolution of the integer nucleon number problem, 
and evidence relating to phase transitions in finite nuclei, are discussed 
from the experimental point of view and interpreted within the frame- 
work of the interacting boson model. 
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It has been known for decades that nuclear properties change, often dramatically, as 
a function of N and Z. However, the possibility of true nuclear phase coexistence and 
phase transitions in the evolution of nuclear structure as a function of nucleon number, in 
the sense of conventional condensed matter systems, has generally been discounted. The 
reason has to do with the finite nature of atomic nuclei and the fact that they contain 
integer numbers of nucleons. In order to discuss the concept of phase transitions one 
needs to identify a control parameter as well as an order parameter. If the fluctuations in 
the order parameter are small, that is, if the data all lie along a compact trajectory as a 
function of the control parameter, the existence of a phase transition would be signaled 
by a critical point where the order parameter is discontinuous. 

The difficulty in applying this approach to nuclei lies in the obvious fact that nuclei 
contain integer numbers of nucleons and, therefore, their properties change discretely as 
a function of N and Z. This is illustrated in fig. la where a typical collective observable 
in even-even nuclei, the energy of the first state, is plotted against neutron number 
for the Sm (Z=62) isotopes. Clearly, an interesting (and well known) change in structure 
is occuring but the abscissa values are, by definition, discrete. Therefore, regardless of 
what the data do, the integer nucleon number requires that, at best, one can only connect 
adjacent points by straight line segments and one can only define differences in properties, 
not derivatives. 

It is the purpose of this Letter to discuss phase transitions in finite nuclei and to 
suggest a resolution of the integer nucleon number problem. When one considers low 
energy nuclear structure, nuclear models offer the flexibility of having one or more con- 
tinuous parameters that can serve as control parameters, allowing the study of critical 
phenomena. The situation in actual nuclei is different. In contrast to phase transitions 
in a speciflc nucleus, where excitation energy (or temperature) can be introduced as a 
control parameter, structure at low energy evolves from nucleus to nucleus as a function 
of and Z, which we have seen are not useful as control parameters. However, we will 
suggest an empirical quantity which is continuous and which does correlate extremely well 
the structural changes with N and Z across large regions of nuclei, and then apply this 
approach to the A = 150 regionj^, ^, |[. Finally, a theoretical analysis with the Interacting 
Boson Model (IBM) |^ shows why sharp transition regions are indeed an expected feature 
of structural evolution. 

We noted in flg. 1 that the data from a single element do not allow a discussion of 
phase transitional behavior. The data for collective nuclei in an entire region (50 < Z < 
66), shown in flg. lb, only exacerbate the problem. The abscissa values remain discrete, 
but now another prerequisite for a phase transition disappears, namely, the absence of 
fluctuations. The scattering of the data obliterates any evidence of sharply discontinuous 
behavior. 

How can we get around this situation in finite nuclei? A possible answer, which we 
propose here, is to choose a qualitatively different quantity to play the role of a control 
parameter, one that is at least potentially continuous. Of course, even if the abscissa 
points become continuously distributed, we need to produce an approach in which the 
fluctuations in the data are small. To do so, consider flg. 2a-c, which should be read as 
successive frames. In this flgure, the same £(45*") data as in flg. 1 are plotted, not against 
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N, Z, or A, but against another collective observable, the energy, E{2i), of the first 2+ 
state. The top left panel is no better for discussing phase transitional behavior than is 
fig. la. The data points are discrete both vertically and horizontally. Although a nucleus 
can, in principle, have any 2f energy, the nuclei (Sm isotopes) in fig. 2a have specific and 
discrete values since their properties change rapidly from isotope to isotope. 

However, when the data for additional elements are added (fig. 2b, c) we see a behavior 
that is qualitatively different from fig. 1. The distribution of points as a function of E(2^) 
successively fills in, yielding in fig. 2c a nearly continuous distribution. Next, we note that 
E(2^) clearly correlates nuclear equilibrium properties extraordinarily well. The data for 
different elements lie along essentially identical traiectohes and thus the ensemble of data 
also lies along a single compact curve, with very small fluctuations. This thereby enables 
a discussion of the trajectory and a potential interpretation in terms of phase transitions. 
Of course, other observables that reflect the equilibrium configuration could also have 
been chosen. 'E{2f) is preferred, however, since it is well and accurately known in a large 
number of nuclei and is easy to measure in newly accessible nuclei. 

Observables such as E{2f), E{4:f), separation energies and other measures of structure 
cannot, of course, rigorously speaking, be considered as control parameters since they are 
not independently variable as is, for example, the temperature in a condensed matter sys- 
tem. Nevertheless, E{2i), de facto, plays a similar role to a control parameter, and fig. 2c 
suggests evidence for true phase transitional behavior as seen by the nearly discontinuous 
change in slope from 2.00 to 3.33 at a specific value of E(2^) which we denote Ec(2|^). 
The expanded view of the rotational region in fig. 2d clearly shows the different slope for 
these nuclei. The change in slope occurs at Ec{2f) ~ 120 keV, although this value may 
differ in regions where other shells are filling. 

Figure 2c itself is not a new result. We have discussed the E(4^) - E{2f) and related 
correlations in ref. and even broached the subject of phase transitional behavior. The 
correlation in fig. 2c has been discussed theoretically^ in the context of the 1 /N expansion 
method for the IBM, and extended^, |T0| empirically to observables for intrinsic excitations 
and to odd A and odd-odd nuclei. What is new here is the explicit discussion of the process 
of reaching fig. 2c as a way of resolving the finite nucleon number problem and identifying 
phase transitional behavior. Of course, since nuclei are finite systems, the phase transition 
is smoothed out over a narrow region of 2^ energies. 

Our next point is to relate this phase transitional behavior to the recent evidence [^, |^ 
for phase coexistence in ^^^Sm (not merely the familiar shape coexistence originating from 
an intruder state mechanism). Theoretical analysis^] of ^^^Sm points to two co-existing 
phases, a deformed ground state band and a spherical anharmonic vibrator spectrum built 
on the 0^ level. Phase coexistence must occur at the critical point and so it is reassuring 
that the 2i energy of ^^^Sm (121 keV) occurs at the same energy as that where the change 
in slope from 2.00 to 3.33 occurs in fig. 2c. 

Thus, two different perspectives - phase coexistence in ^^^Sm and the relation of yrast 
energies across the region of nuclei - give evidence for a phase transition from spherical 
to deformed structures, near A = 150. Note that, while phase coexistence occurs in a 
specific nucleus (^^^Sm), the phase transition does not characterize a single nucleus, or 
the isotopes of an element, but is a property of, and only definable in terms of, an entire 
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region. We stress here that this type of phase transition and phase coexistence is different 
from the shape coexistence picture known in other regions |[TT| , |12]. Here the structural 
changes develop within the context of a single shell and do not involve an intruder state 
mechanism. 

Other observables reflect these rapid structural changes. Two nucleon separation 
energies, S2n, are sensitive to changes in equilibrium configurations. In fig. 3a we show 
the empirical values of for the 50 < Z < 66 nuclei (for all collective nuclei - those with 
R4/2 > 2.05). The values have a well known, essentially parallel, shift in values for 
each successive Z. To compare values for different elements on the same scale, we therefore 
shift the separation energies for each Z by a constant amount chosen to give equal S2n 
values at N = 88 for the N > 82 shell and at N = 76 for the N < 82 shell. 

For N>82, the results for S2n in fig. 3 striking as for the Af energies in fig. 

2c. The behavior is compact, with small fluctuations and a sharp break in trajectory. 
This break occurs at a slightly lower 2f energy than the slope change for £(45*") in fig. 
2c. Apparently, the S2ri values display a different dependence on the overall shape of the 
nuclear potential than low spin yrast levels (which presumably are most sensitive to the 
details of the potential near its minimum). For the lighter mass shell, N<82, there are 
greater fluctuations and a gradual structural change but no evidence for a sharp phase 
transition. 

In fig. 3b we show B(E2: 2^ — > Oi) values as a function of 5'2„. The sudden break of 
trend at a specific value of S2n illustrates the point that £(4^'') is not a unique measure 
of the structural transition: B{E2) values provide additional evidence for it. 

The question arises as to why nuclei should behave in this way. As we have noted, it 
is easier to look at phase transitional behavior in a model since the model parameters are 
inherently continuous. Here, we consider the IBM since the analysis is simple but similar 
results characterize the Geometric Collective Model (GCM) |15|. We use the Extended 



CQF Hamiltonian[16, 17 



H = end - kQ ■ Q, (1) 

where Q = {s^d + d^s) + x{d^dY, which has parameters e, k, x and the boson number 
Nb- Performing a systematic analysis of parameter space, IBM calculations of E{4f) 
versus -^(2^^) reproduce the results of fig. 2c: they follow a slope of 2.00 above Ec{2f) for 
virtually any choice of e, x, and that gives R4/2 = 2.05 - 3.15 as long as k, is constant 
p!8| . Indeed, the value of k determines the intercept €4 [and Ec{2f) as well]. Were the 



data different from fig. 2c (e.g., scattered, or following a slope other than 2.00) the only 
way this behavior could be reproduced would be to separately adjust k for each nucleus, 
a clearly unlikely scenario, and one that is inconsistent with microscopic analyses of IBM 
parameters [|19] . 



We now study the behavior of eq. (1) through the intrinsic state formalism [pO|, |21 
computing the energy surfaces corresponding to different parameter values. It is conve- 
nient, for this purpose, to re-write eq. (1) with a control parameter ^ = (1 + e//t)^^. The 
resulting scaled Hamiltonian has the form 
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(2) 



For X = "V^ /2, < < 1 maps the transition from U(5) to SU(3). 

In fig. 4a we show the IBM energy surface corresponding to the classical limit of 
eq. (2), for = 10 and x = -V^f^, against ^. The figure immediately shows the key 
point, that the location of the minimum in the energy, Pmin, changes suddenly, for these 
parameters, at a particular ^ value, from Pmin = to a large finite value, as indicated 
by the dark line cutting through the contour plot. There are virtually no ^ values for 
which intermediate /3mm values result. This is consistent with the specific calculations in 
ref. 0]. The IBM indicates that nuclei change abruptly from near spherical to deformed 
at a critical value ^crit- The evolution of the energy surface can be seen as the competi- 
tion between two minima, spherical and deformed, rather than a gradual evolution from 
spherical to weakly deformed to large deformation. Such a level crossing scenario is, in 
fact, characteristic of a first order phase transition. 

The qualitative behavior of these results is not sensitive to boson number nor to 
variations in x- We show this in fig. 4b, which gives the values of the location of the 
lowest minimum of the energy surface, as a function of ^ for a set of and x values. The 
curves all show the same behavior independent of the parameters: /3mm is zero for small ^ 
values, and then rises extremely rapidly to a saturation value within a very narrow range 
of ^ values (which define a ^crit for each Nb and x)- 

In fig. 4c, we show energy surfaces as a function of P for IBM parameters applicable 
to each of the Sm isotopes. These surfaces range from near vibrator shapes for i46,i48gj^ 
to softer in ^^°Sm, to the unique coexistence nucleus ^^^Sm where two shallow minima 
occur, to i54,i56gj^ where the prolate minimum is much deeper, thus localizing these latter 
nuclei at large positive (3 values. The actual minimum in the energy surface occurs only 
for /3mm ~ or large positive /3mm. 

In summary, we have shown that finite nuclei can, and do, exhibit true phase transi- 
tional character by proposing a way to resolve the integer nucleon number problem. We 
have shown that one can identify an empirical quantity, E(2^), that is nearly continuous 
and in terms of which other quantities, such as E(4^) or S2n, follow simple, compact tra- 
jectories with small fluctuations even for large regions of nuclei. We have seen that these 
observables have distinct anomalies at similar values of E(2^) at which phase coexistence 
in ^^^Sm has been recently identified. Other observables such as may also be used to 
correlate structural changes. Through a model, such as the IBM (the GCM gives similar 
results), we have associated the order parameter with a physical quantity, the deformation 
Pmin, at which the potential has a minimum. In a spherical-deformed transition region 
in the IBM, /3mm takes on two characteristic values, zero and then an abrupt change to 
a near saturation deformation. It appears that nuclear structural evolution in this mass 
region entails two basic phases (spherical and deformed) rather than a gradual softening 
(with valence nucleon number) traditionally associated with the onset of deformation in 
nuclei. Although this view is unconventional, the analysis of the IBM suggests that it is 
nearly unavoidable and a basic feature of structural evolution. 
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Figure Captions 



1: E{Ai) against neutron number for collective nuclei, that is, nuclei for which R4/2 = 
E(4^-)/E(2^-) > 2.05. a) Sm, b) the 50 < Z < 66 region. Data from ref. [2] 

2: E{4f) plotted against E(2l). Panels a-c), show the same data as in fig. 1 for 
sequentially more elements. Panel d) is an expanded view of the rotational region, 
with all the data points from panel c) with E{2f) < Ec{2f). 

3: (a) Separation energies as a function of E{2f). data from ref. [13]. (b) 
B{E2 : 2t Of) values as a function of S2n- B(E2) data from ref. [14]. 

4: Classical limit analysis of IBM calculations that reproduce fig. 2c. a) Energy surface 
as a function of ^ [see eq. (2)] near the critical point for = 10 and x — -V7/2; 
b) location of the minima, Pmin, as a function of ^, for several values of and x; 
and c) The IBM energy E(/3) for the C, values and boson numbers corresponding to 
the Sm isotopes (sec ref. [22]). The minima occur only for /3 = or large finite 
values. There is no gradual evolution of P from to saturation levels. 
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